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1 Introduction 

This paper develops the Wilson line approach to high energy scattering proposed in [1] to 
the case of quadrupole and double dipole operators in the next to leading order (NLO). 
Such operators naturally appear when one studies amplitudes for diffractive processes with 
the production of 3 or 4 particles in the Regge limit. Moreover, the quadrupole operator 
enters the definition of the Weizsacker-Williams gluon distribution [2], [3], [4] which gives 
the Foc.k space number density of gluons inside dense hadrons in light-cone gauge. One can 
find the NLO evolution equation for the operator necessary for the Weizsacker-Williams 
gluon distribution differentiating the quadrupole equation obtained in this paper. This 
result is going to be presented in a future work. 

In the Wilson line approach to high energy scattering [1] the amplitudes are convolu¬ 
tions of impact factors and a Green function. The impact factors describe the decomposi¬ 
tion of the colliding particles into quarks and gluons while the Green function is responsible 
for the interaction of these quarks and gluons with the quarks and gluons from the other 
colliding particle. In this framework such fast-moving partons are depicted as Wilson lines 
with the path going along their trajectories. Hence, the corresponding Green functions are 
the operators constructed of the Wilson lines. These operators obey the evolution equa¬ 
tions with respect to the rapidity divide. This rapidity divide separates the gluon field into 
the fast quantum one and the slow external field of the other particle, through which the 
current quark or gluon is propagating. 

In the most thoroughly studied case of a virtual photon splitting into quark antiquark 
pair, the corresponding Wilson line operator is a color dipole. The evolution equation 
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for this operator is known as the Balitsky - Kovchegov (BK) equation [1], [5]. The NLO 
corrections to this equation were calculated in [6], [7], [8], [9]. Another interesting case is 
application of this formalism to a proton. The proton has baryon color structure and can 
be described as a 3-quark Wilson loop operator (3QWL). The evolution equation for this 
operator was calculated in the leading order (LO) in [10] and in the NLO in [11]. The latter 
calculation was based on the NLO hierarchy of the evolution equations for the Wilson lines 
with open indices [12] and the connected contribution to the 3QWL kernel [13]. These 
results were also obtained in the Jalilian-Marian, Iancu, McLerran, Weigert, Leonidov and 
Kovner (JIMWLK) formalism [14], The hamiltonian equivalent to the NLO hierarchy was 
obtained in [15] and the evolution equation for the 3QWL in [16]. The NLO kernel for the 
evolution Wilson line operators was also constructed in [17]. 

The quadrupole and the double dipole are 4-particle colorless operators. Their LO 
linear evolution equations were derived in [4], [9], [18], [19]. Here the results of [12] and 
[13] are used to construct the NLO evolution equations for these operators and the results 
of [9] and [11] are used to check these equations. 

The paper is organized as follows. The next section contains the definitions and neces¬ 
sary results. Sections 3 and 4 present the NLO evolution equations for the quadrupole and 
the double dipole operators in the standard and quasi-conformal forms. Section 5 discusses 
different checks of the results. Section 6 concludes the paper. 

2 Definitions and building blocks 

The light cone vectors n\ and ri 2 are defined as 

ni = (1,0,0,1), n 2 = ^ (1,0,0,-1), nf = = nin 2 = 1 (2.1) 

and any vector p can be decomposed as 

p + = p- = pri2 = - (p° + p 3 ) , p+ = p~ = pn i = p° — p 3 , (2.2) 

p = p + n\ + p~ri 2 + p±, p 2 = 2p + p~—p 2 , (2.3) 

pk = p 9 k^ L = p + k~ + p~k + — pk = p+k- + p~k + — pk. (2.4) 

For brevity the following notation for traces is used 

triUiUl.-.Uklh t ) = U ijt ... Htl (2-5) 

where 

Ut = U (n , V ) = Pe i9 ■/■-» b * r+ ’ rldr+ , (2.6) 

and b~ is the external shock wave field built from only slow gluons 

K = f 7WT4 e "^" (p) W - P + ). (2.7) 
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The parameter i] separates the slow gluons entering the Wilson lines from the fast ones in 
the impact factors. The field 

6^ (r) = 6 _ (r + , r)n 2 = 5(r + )b (r) n 2 . (2.8) 


The coordinates ?q, 2 , 3,4 denote the quarks, and f*o,r 5 are the coordinates of the gluons. In 
intermediate formulas the coordinates ffjj will also be used. The SU ( N c ) identities 

Ut a = 2 tr(t b U 4 t a U\), 0 nini = 2 S l A - (2.9) 

are necessary to rewrite the SU(N C ) operators only through the Wilson lines in the fun¬ 
damental representation. For a generic operator O the rapidity evolution equation has the 
form 

— (0) = (Klo®0) + (K nlo ®0). (2-10) 

where Klo ~ and Knlo ~ The (...) brackets were explicitly written to denote that 
the calculation was performed in the shockwave background. Hereafter they will be often 
omitted to avoid overloading the notation. The BK equation in this notation reads [1] 

—r r- = Trk f dro AAA-2 ( U 2toUot! - iV c U 2tl ), (2.11) 

dt] 2i r- J r 10 z r 20 z 

where r t j = Vi — fj. The LO quadrupole evolution equation reads [4] 

= 4 ^ j ^ 0 { n 0 V 40 2 ^ Ul0tU ° 2t34t + U 4t0 U l2t30t - (0 1 = 0 4)) 

+ - ^ ^ 2 (Ul0tU 0 2t34t + U 2 toUi 0 t34t — (0 —> 1 = 0 —> 2)) 

rio r 2 o 

r 2 4 2 

— 2 - 9 (UiotU 02 t 3 4t + U 30 tU 04 t 12 t — (0 —>■ 4 = 0 —> 2)) 

2r 20 ? 40 

_ 9 -> 2 -> 2 (U 4 toUi2t 3 ot + U 2 toU 34 tiot — (0 —)■ 1 = 0 — > 3)) + (1 <-)■ 3,2 <-)■ 4)1 . 
2^10^30^ J 

( 2 . 12 ) 


Here (0 —> 1 = 0 —> 4) stands for the substitution U$ —> U\ or Uq —> U 4 , which gives the 
same result. In addition (1 f> 3,2 o 4) means that one has to change f\ «-)• f, 3 , r 2 <->• r 4 
and U\ <-)• U 3 , U 2 «->■ U 4 . We will also need the LO evolution equations for the double 
dipole, sextupole and the dipole-quadrupole product. All these equations follow from the 
LO hierarchy [1] directly. 


dU 12t U 34t a s f 
Or, = 4^ J dn 


r 13 


r 23 


ru 


Tio 2 r 30 2 


T20V30 2 


n 0 2 r 40 2 


+ 


r-24 


T 2 o 2 r 4 o 2 


x (U 2 ti4t3 + U 2t34 ti - U 2 t 10 t3 4 to — U 2 t 04 t30tl) + u 4t3 + U 2 t 1 4 3 


d ij 


dr/ 

(2.13) 
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“ — -|^2 j {(Ugt76t02t34t 1 + U 0 t 12 t34t06t7 — (0 —> 7 = 0 —> 6)) 

f r 16 2 _ fj 7 2 \ / f 2 7 2 _ ^* 26 2 \ 

X \r 0 i 2 r 0 6 2 r 0 iM) 7 2 y \ro 2 2 ro 7 2 r 02 2 r 06 2 J 

x (Uot76t02t34tl + U 0 t 3 4t!2t06t7 - (0 —> 7 = 0 —> 6 )) 

+ d » 3,2 « 4)} + U 76 , + U 12t31t ^. (2.14) 


_____ ^ ro {___{jJ of3 4t 5 U 2 t ofi ti + Uot56llU 2 l34to — (0 —45 = 0 —>2)) 


dij 


>02 2 n)5 2 


47T 2 t 
r 15 * 

— ^ 2 -> 9 (Uot56tlU 2 t 3 4to + U 6 t 0 U 0 t 12 t34t5 — (0 — > 5 = 0 —> 1)) 
?01 H)5 

foe: 2 

_ - 2- 2 (Uot34t5U 2 t06tl + U 0 t 1 U 2 t 3 4t 56 to — (0 —> 2 = 0 —> 6)) 

?’02 r 0 6 

f ,2 

+ - 2 -> 2 (U 6 t 0 U 0 t 12 t34t5 + U 0 t 1 U 2 t 3 4t56to — (0 -> 1 = 0 -> 6 )) 
? 01 H)6 

—* o 

H2^ 


+ 


(U 2 toU 0 t 3 4t56tl + U 0 f4U 2 t34156+0 — (0 —> 1 = 0 —> 2)) 
+(1—>3—>5—>1,2—>4—>6—>2) + (l—>5—>3—>1,2—>6—>4 


rm 2 r 0 2 2 


■ 2 )}. 

(2.15) 


Here 1 —> 3 —> 5 —> 1 stands for permutation, i.e. one has to change f\ —> f 3 , 7% —> r* 5 , 

f 5 -> r=i and Hi -> U :i , H 3 -> H 5 , H 5 -> Hi. 

For the self and the pairwise NLO interactions one can take the results of [12] while 
the triple-interaction diagrams were already calculated in [13]. The results of these papers 
were derived using sharp cutoff on the rapidity variable. Since this paper is devoted to color 
singlet operators one can drop the kernels which vanish acting on the colorless operators, 
as was shown in [15]. The rest reads 


d (Hi) 


di] 


NLO 


= df 5 df 0 Jnio 5 [if ad ' e \{t d t e }U 1 ty i -if ad *(t a U 1 {t*t e '})l] 


Tjdd 1 (rree 1 _ Tjee\ , a s N c f dj^ 

x U 0 [U 5 U 0 ) + 47r 3 J ^ 


(Uf - Hf )(HHiH)^ln ( ) , (2.16) 


T a b\/ + a TT_ +b\j 


15 


P = 


11 2 rif 

3 ~~ 3 Nr 


ft In ^ = 


11 2 n f 

T _ 3 Wr 


In 


( d 2 

V 4e 2 ^’( 1 ) 




67 7T 2 10 iif 

+ y _ y _ y ah 


(2.17) 


n/ is the number of the quark flavours, fi 2 is the renormalization scale in the MS-scheme 
and Jijkim = J{n , fj,r k ,f t , f m ) 


«7l2305 = 


H)H52 , 2 (r 0 ir 03 ) (ro 5 r 25 ) 


roi 2 r 05 2 T 25 2 


+ 


Tor 2 7U3 2 To5 2 f 2 5 2 


x In 


r 03 

y 2 


2 (roir 03 ) (yy) 2(roiro 5 ) (r 25 r 35 ) \ 

r 0 i 2 y 2 y 2 y 2 y 2 y 2 y 2 y 2 y 

(2.18) 
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This function has the properties 


T _ T T (j 5 iroi) , (m 2 

^ijk 05 JjikbOi ^11105 -» 9^1 -* 9 

roi z ro5 z ris z \ r i5 


(2.19) 


0{Ui)UU 2 ) l k = a^_ r df ^ ^ +A2 + As) + ^_ f ^ (Bi + . (2.20) 


=[(t a C/i)|(C/ 2 t 6 )L + (t a t/ 2 )i(C/i^] - C/ 0 ee ')4L 12 

+ ^nfL\ 2 tr{t a U^t b {ul — C/|))l , 


( 2 . 21 ) 


where L,j = L(r),r) ) and L c j- = L g (ri,fj) were introduced in this form in [11] 


L 12 — 


roi 2 ^ 2 5 2 - rb 2 2 ri5 2 


'Tbi 2 r25 2 + r 0 2 2 ri 5 2 ) + r 05 2 


n 2 2 r 02 2 ri 5 2 + r 0 i 2 r 25 2 


n 2 2 / 1 _ 1 X\ ] / r 0 i 2 r 25 2 \ 1 

8 ro 5 2 \ro 2 2 ri 5 2 T 0 i 2 r 2 5 2 /. Vtls 2 ^ 2 / 2 t ? 0 5 4 ’ 

T q _ 1 f r 0 2 2 ri 5 2 + r 0 i 2 r 25 2 - r 05 2 ri 2 2 / r 02 2 ri 5 2 \ 

12 r 05 4 \ 2(r 02 2 ri 5 2 - r 0 i 2 r 25 2 ) V^oi 2 ^ 2 / 


( 2 . 22 ) 


(2.23) 


These functions have the unintegrable singularity at ro 5 = 0, which is canceled by the 
subtraction in the color structure. They are symmetric conformally invariant functions 

r(g) _ r(?) _ r (?) I 

— j* — i? ho-H>r 5 - 

•4 2 =4(t / 0 - C/i)^(t / 5 - ^ 2 ) ee ' ^i[f ade '(t d Uit a y i (t e U2) l k - f ade (t a U 1 t d ')i(U 2 t e '){}J 12 i 05 
+ i[f ad ' e '- / ade (C/it d 'H(t a C/ 2 t e ')i] J 12205 } , (2.24) 

Az =2 Uf {i[/ ad ' e Vit a )K^ e ?7 2 )i - f ade (t a U\)^(U 2 t e 't d 'Yfcl (C /5 - C/ 2 ) ee ' 

X (</ 22 105 + </ 2 1205 — </l 22 05) + (<Cl2105 + <^11205 “ <^21105) 

x i[f ade '(t d t e Ui) :, i (U 2 t a ) l k - f ade (Uit e 't d ')i(t a U 2 ) l k \(U 5 - C/i) ee '} . (2.25) 




X |(C / 5 - Ui) ab i[f bde (t a Uit d )i(U 2 t e y k + f ade (t e Uit b )i{t d U 2 ) L k ] (_ _L 


15'25 '15 


+ (u 5 - c/ 2 ) a6 z[/ 6de (c/it e )j(/ a c/ 2 t d )i + r rfe (/ d c/i)j(t e c/ 2 t 6 )y 

£2 =/3 (2C/ 5 - [Cl - C/ 2 ) ab [(C a C/i)KC/ 2 t 6 )i + ([/i^(C a [/ 2 )i] 


,6W i Z' (rL5r 25 ) 1 


—* 2 ~* 9 —* 2 

' 15' 25 ' 2 5 


f>2 zf 2 

'15'25 


(n-Ai) ln (fut + J_ ln ('M') + J_ to CfS 

l r.2 ) ^ 2^2 V ^2 1 ^ 0 , 1-2 l ^*2 
\ / z, 15 \ ' 12 / z/ 25 \ ' 12 


(2.26) 


(2.27) 
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m: 


la 


NLO 


= / dr 5 dr 0 


X {f cde ((t a Ui)i (t b U 2 y k (U 3 t c ) n m (Up - Th) M (U b - u 2 ) 


ad( 


\be 


- (u.ni mx ( t c u 3 r m (u 0 - ux a (u 5 - c/ 2 ) e Vi2305 

+ f ade [(Uit a )i (t b U 2 ){ (t c U 3 )l (U 0 - U 3 ) cd (U 5 - U 2 ) be 

- (fUXi (U 2 t% (U 3 t c ) n m (Uo - U 3 ) dc (U 5 - C/ 2 ) e V32105 
+ f bde [(t a Ui)i (U 2 t% (t c u 3 r m (Uo - c/i) arf (c/ 5 - U 3 y e 

- (u.ni (t b u 2 y k ( u 3 n n m (Uo - uy da (u 5 - u 3 r]Ji32 05 }. 

(2.28) 


We will also need the following functions. The function M 3 - k = M (r), rj,f k ) was introduced 
in [ 11 ] 


M 2 3 — ~(Jl2205 + J23205 ~ <^ 13 2 05 “ «/ 222 05 ) 


1 


^ 2z? 2 ^ 2 ,v 2 z? 2zt 2 Zf 2 

r \2 r 23 ri 5 ? 2 3 H )3 r \2 r 13 

I — o 


In 


r 02 
V 2 U 


4r 0 i 2 r 35 2 V ro2 2 r 2 5 2 r 0 5 2 r 2 5 2 ro 2 Vo 5 2 r 0 5 

It was also introduced as M 2 in [20]. It has the property 

= •'//!'• 

The functions = L(fi,fj) and = M(fi,fj) were introduced in [11] as well 

1 1 


(2.29) 


(2.30) 


L 12 = -(Ml 2 - M 2 n ) = 

2 V 1 2 ; 8r 05 2 


m 12 = 1(m 2 11 + mP) = i || ? 


i? 2 ^ 2 

n 2 ^05 


nji 2 To 2 2 r 1 5 2 f 25 2 r 0 i 2 r 2 5 2 r 02 2 fi 5 2 


In 


2 ^ 2 
Hu r 23 


ri5 2 ro2 2 


zf 2zt 2 

ri 2 r 0 5 


1 


1 


'doi 2 To 2 2 ri5 2 r 2 5 2 r 0 i 2 r 25 2 r 02 2 r\r 2 


In 


Zf 2zf 2 
rp\ r q 2 

ri 5 2 r 25 2 


(2.31) 


(2.32) 

Here Lij is conformally invariant. Moreover, L^j is antisymmetric w.r.t. both 5 -H- 0 and 
i -H- j transformations while My is antisymmetric w.r.t. 5 -B- 0. One can also combine all 
the terms ~ f5 into = Md(fi,fj) 


Mf 2 = M 
12 2 



+ 


^ 2-2 

^ 2 ln 

r 12 





(2.33) 


The NLO BK kernel reads [9] 

(K NLO 9 U 12t ) [dr 0 J^- 2 { Mf 2 - N C In f In f §£) 

4vr-V ri 0 -r 20 2 ( V r io / V^o / 

a 2 r 

x (U 2 toU 0 ti - iV c U 2 t i) + ^4 J drodr 5 {L 12 (\J 0 ^V 2io U 5n -(0« 5)) 

+Li 2 ((U 0 t 52 t 0 5 ti - U otl U 2 t 5 U 5 to — (0 —> 5)) + ( 0 f> 5)) 
-2n f L\ 2 (tr(t a Uit b ul)tr(t a U 3 t b (ul - u\)) |(5h 0))}. (2.34) 
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3 Construction of the kernel 


First, one has to discuss the singularities of the building blocks from the previous section. 
All the ultraviolet (UV) singularities in (2.16), (2.20), and (2.28) were removed by the 
renormalization. It means that these expressions converge at r*o = ri, 2 , 3 , 4,5 and F 5 = 
A),1 , 2 , 3 , 4 - hr particular, the functions J in A 2 (2.24), A 3 (2.25), and (2.28) are convergent 
at these points, which ensures UV-safety of these expressions. However, the function J 11105 
in the first line of (2.16), has the UV singularity at ro = F 5 = r\. As in (2.22) this singularity 
is removed by the subtraction in the color structure. 

Nevertheless, these expressions have infrared (IR) singularities, which appear as both 
f*o ,5 — > 00 . Indeed, changing the variables e.g. as ?’o = ut, r$ = ut, t = 1 — t, one faces a 
logarithmic singularity integrating w.r.t. u first 

/ dr 5 dr 0 J 12305 = / #5#0 dt£ du + In (f) +' O (4)) . 

(3.1) 

Hence this double integral is ill-defined and requires either regularization or definition in 
terms of the iterated integrals. To understand how to correctly treat the IR singularities one 
can either return to the diagrams and keep the regularization, or calculate the known dipole 
equation and fix the definition from there. The latter way is attempted here. Assembling 
BK kernel (2.34) from (2.16-2.20), one can see that all the /3-functional terms go to m£, A\ 
( 2 . 21 ) reshapes to the terms Li 2 ,L\ 2 , the Wilson line operators from (2.16), (2.24-2.25) 
depending on both r*, and fo give the term ~ L 12 after the symmetrization 


A(U)F(r)-> [AF} sym 

[4 + 4 (0 H 5)] [F + F (0 H 5)] + [A — A (0 ** 5)] [F - F (0 o 5)] , o 
- , • ( 4 - 2 ) 


Next, B\ (2.26) gives one half of the double logarithm contribution. All the remaining 
terms are to be equal to the other half of the double logarithm contribution. They read 


A K = 


ot 2 s N c 

167T 4 


J dr 0 (ir5{(U2t5U 5 t 1 - U 0 t 1 U 2 t 0 ) (T22105 + ^11205) 


+ (<321105 — T12105 + <3l2205 — <321205) (2V c U 2 ti — UotiUjjto — U 2 t5U 5 t!)}. (3.3) 


This term is IR safe. The second line is the product of 2 expressions symmetric w.r.t. 
0 n 5 permutation. Therefore one can set U 5 —> Uq there. In the first line there is a 
product of 2 expressions antisymmetric w.r.t. 0 <-)• 5 permutation. Hence, one could add 
and subtract V c U 2 t 1 in the first brackets and write 

AA — > J drodr5{(N c U 2 n — U 2 t 0 U 0 ti) (<322105 + <3mo5) 

+ (<321105 + <3l2205 — <3l2105 — < 321205 ) (3V c U 2 t 1 — U 0 t 1 U 2 t 0 )} (3.4) 

= - ^/ dro J dr 5 (N c U 2tl - U 0tl U 2t0 ) 

x{ J 22105 — <321205 + <311205 ~ <3l2105 + <321105 + <3l2205)< (3.5) 
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One could understand the latter integral as an iterated one. Then, using the integrals 


( Jijk05 JikjOh) 


roiroj 

^2zf2 

Oi Oj 


T0i r 0k 

i?2j?2 

Oi Ok 




j df 5 Jijj 05 = 0, (3.6) 


one could get the other half of the double logarithm term in the BK kernel 


A K 


ot 2 s N c 

87T 3 


dro (iV c U 2 ti — U 0 tiU 2 t 0 ) 


r 12 


ri 0 2 r 20 2 


In 




(3.7) 


Although such treatment gives the correct result, it does not take into account the IR 
singularity of J. Indeed if one introduces the dimensional regularization into (3.3) then one 
gets 

A7i d d r 5 (N c U 2tl - U 0tl U 2t0 ) 

x (<722105 — <721205 + <7ll205 — <7l2105 + <721105 + <7l2205)< (3.8) 


However, in the dimensional regularization the integral f d d r§Jijj 05 would be ~ e rather 
than 0 and the double integral 


J d d r 0 d d r 5 J ijj05 = 2vr 2 C(3) 


(3.9) 


because the second integral w.r.t. vq has an IR divergence as ro —> oo and starts from 
7. Therefore if one wants to integrate the coefficient of U 0 t 1 U 2 t 0 w.r.t. d d rs, one has to 
keep the result in the dimension d without expanding the series. At the same time the 
coefficient of IV c U 2 ti gets the doubled contribution since 


df 0 


r 12 


no 2 r 20 2 


In 


n 2 
Tio 2 



4vr 2 C(3). 


(3.10) 


Thus, the result depends on the regularization. Such an ambiguity is the consequence of the 
fact that the initial expressions do not have the IR regularization. To avoid this ambiguity 
one needs the evolution equations for Wilson lines (2.24-2.25) with the IR regularization. 
Alternatively, one can write them in the form where the terms which do not depend on 
both U£ b and t/g 6 are integrated w.r.t. the coordinate of the other gluon. 

In this paper the procedure discussed in (3.2-3.7) is used. Technically it means that 
for the terms ~ U£ b UQ b ', the gluons are treated equally and the kernel is represented in 
the form of symmetrized sum (3.2). In the terms depending only on U£ b or only on t/g fe , 
the integration order is fixed as f dro f dr§ or f dr 5 f dro correspondingly and the integrals 
are understood as iterated. As a result, one can take the inner integral via (3.6). The 
terms independent of U£ b and I/g b are also symmetrized according to (3.2) and in them 
the substitution Jijj 05 —> JjjiOb is made. This substitution can be understood as follows. 
First one drops the terms with f dr $ Jijjo^ ■ They vanish (3.6) if one treats the integrals as 
iterated. Next, one adds the totally antisymmetric w.r.t. (5 -B- 0) terms Jjjiob- These terms 
vanish if they are integrated w.r.t. r*g and F 5 in the double integral. After that the first 



integral in (3.6) is enough to calculate all the integrals. Again, I stress that although such 
treatment gives the correct dipole result (as well as the evolution equation for the baryon 
operator coinciding with [ 11 ]) it involves the cavalier treatment of the IR singularities. 

Taking the contributions of the self-interaction of one Wilson line (2.16), the con¬ 
nected contributions of 2 Wilson lines (2.20) and the connected contributions of 3 Wilson 
lines (2.28) with the appropriate charge conjugation, and using the integration proce¬ 
dure described above, one can write the full NLO evolution of the quadrupole operator 
tr {XhUjlhUj) = U 12 t 34 t as 


o? f o? f 

(K nlo ® U 12 t 34 t) = ^J dr 0 df 5 (G s +G a )+^ j df 0 (Gp + G), (3.11) 

Following (3.2) the 2-gluon contribution can be decomposed into the product G s of the 
symmetric coordinate and color structures and the product G a of the antisymmetric ones 
w.r.t. 0 5 transposition, i.e. the substitution fo -H- rb and Uq <->• C/ 5 . After color 

convolution and integration w.r.t. vq or 7/5 of the contributions which do not depend on 
the other variable one comes to the 1-gluon part. It contains the contribution proportional 
to /3-function Gp (/3 = ^ — |^f) and the rest G. One can see that all the G’s separately 
vanish without the shockwave, i.e. if all the U —> 0. 

Doing the same for the double dipole operator tr (UiU 2 ^)tr(U 3 U^) = U 12 tU 34 t, one 
can write its full NLO evolution equation as 

{Knlo <8> U 12 tU 34 t) = U 12 t(Al/VLO <8> U 34 t) + G za \{I\nlo <8> U 12 t) 

+ J drodrs (G s +G a )+^Aj J dr ,b (G^ + G). (3-12) 

Here the NLO dipole kernel is written in our notation in (2.34), G s (G a ) is the product 
of the coordinate and color structures (anti)symmetric w.r.t. 0 O 5 transposition, Gp is 
proportional to /3-function and G is the remaining contribution with 1 gluon crossing the 
shockwave. 


3.1 Quadrupole 

We start from the product of the symmetric structures 

G s = G s i + rifGq + G s 2- (3.13) 

G ,1 = ({U 0 t34ti 5 t02t5 - U 5 t 0 U 2 t 5 U 0 t 34 tl - (5 -> 0 )} + (5o 0 )) (L 12 + L 32 - L 13 ) 

+ ({U 0 tl 5 t 02 t 34 t 5 - U 0 t 5 U 5 t 1 U 2 t 34 to — (5 “^ 0 )} + (5 <R- 0 )) + Lu ~ A 42 ) 

+(1^3,2o 4), (3.14) 
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where L was introduced in (2.22). It is a conformally invariant contribution. 


G, ^ ( Pjg«UW _ yafefll _ U 2 , 5 U 0t34tl - (5 - 0 )} + (5 « 0 )) 

x 2 (^12 + r 32 - r 43 ) + 2(^12 + r 44 — l 42 ) 

^/rUotl2t34t5 + U 2 t34tl5t0 U 0 t 5 U 2 t 34 ti , . nU . ( K . . n>\\ 

x(l--- U 5tlU 2 t 3 4tO - (5 -> 0 )} + (5 0 )) 

+(1 h3,2h4), (3.15) 

Here L q is defined in (2.23). 

G S 2=- (U 0 t 34 t52to5ti ~ U 0 t 1 U 2 t 5 U 4 t 05 t 3 + (5f> 0)) (Mf 1 — Mf 4 + Mf — M 2 3 + (5f> 0)) 

+ 9 (Uot 35 t 02 t 54 tl — U 0 t 3 U 2 t 5 U 4 t 15 to + (5 o 0)) (M ^ 4 — Mf 4 + Mf 1 — Mf 1 + (5f> 0)) 

+ 9 (Uoti5to2t34t5 — U 0 t 5 U 5 t 1 U 2 t 34 to + (5 o 0)) (Mf 4 + Mf 2 + (5f> 0)) 

+ 9 (U 0 t34ti5t02t5 “ U 5 t 0 U 2 t 5 U 0 t 34 ti + (5 -H- 0)) (Mf 3 + Mf 1 + (5 a 0)) 

+(1 f>3,2f> 4), (3.16) 

where Mj k is defined in (2.29). Using property (2.30) one can show that G s2 vanishes 
without the shockwave, i.e. when all the U —> 1. Indeed, it is clear from the representation 

2G<j2 = (U 2 t 5 U 5 t 0 U 0t34 ti - U 4 t 0 U 5 tl U 0 t 52t3 + U 0 t 15 t04t52t3 _ U 0 t 34 ti5t02t5 + (5 0)) 

x(M 32 — Mp) 

+ (U 4 t 0 U 5 t 1 U 0 t5 2 t 3 -U 2 toU 5 t 1 U 0 t 34 t5-U 0 t 15 t 04 t52t3 + U 0 ti5t3 4 to2t5 + (5 ^ 0))(M 4 2 —M 2 ) 
+ (U 0 tiU 2 t 5 U 4 t 0 5t 3 -U 4 t5U 5 toU 0 t 12 t 3 + U 0 t 12 t 3 5to4t5 _ U 0 ti5t3 4 to2t5 + (5 ^ 0 ))(M 43 —M 34 ) 
+ (Uot3U 2 t 5 U 4 t 1 5to-U 2 t 0 U 5 t 1 U 0 t 34 t5 + U 0 t 15 t 34 to2t5-U 0 t52t04tl5t3 + (5 ^ 0))(M| 4 — Mf 3 ) 
+ (Uot5U 5 t 1 U 2 t 34 t 0 -U 0 t 3 U 2 t5U 4 t 1 5t 0 -U 0 t 15 t 02 t 34 t5+U 0 t35to2t54tl + (5 <H- 0))(M 2 1 —Mf ) 
+ (Uot5U 5 t 3 U 2 t 04 t 1 —U 2 t 0 U 5 t 1 U 0 t 34 t 5 +U 0 t 15 t 34 t02t5 —U 0 t5 4 ti2t05t3 + (5 ^ 0))(Af| 3 —M| 4 ) 

+ (le3,2f>4). (3.17) 

The contribution which is the product of the antisymmetric w.r.t. 5 <->• 0 parts reads 

G a = G a i + G a2 + G a3 . (3.18) 

G a i=- (U Qt xU 2 tgU 4 t 0 5 t 3 + U 0 t 34 t 52 t 0 5 ti — (5 -H- 0)) (M 2 31 — M | 4 — Mf 2 + M 43 - (5 0)) 

+ 2 (Uot3U 2 t5U 4 t 15 to + U 0 t 35 t 0 2t54ti - (5 o 0)) (M 43 — Mf 4 - Mf 2 + Mf 1 - (5 0)) 

+(1 o 3,2 o 4). (3.19) 
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G a 2 — — (Uot34ti5t02t5 ~ (5 -H- 0)) (-£13 + 2M‘2\ — 2M23 — A/ 23 + M 21 - (5f) 0)) 

+ — (Uoti 5 to 2 t 34 t 5 — (5 -H- 0)) (L 42 — 2M\2 + 2M \4 + M 2 4 — M\ 2 - (5-b 0)) 
+(lo3,2o4). (3.20) 

Here the functions L and Mjj are defined in (2.31) and (2.32). 

G a 3 =— (UgtgUgtiU 2 t 34 to -(5 -h- 0)) (L ±2 + T 14 — 2 M 24 + Af ^ 4 + Af | 2 — (5 o 0)) 

+ - (U5toU 2 t5U 0 t34ti — (5 o 0)) (L 21 + L 23 — 2Afi3 + Aff 3 + Afg 1 — (5 o 0)) 
+(lo3,2o4). (3.21) 

From (2.31) and (2.32) one can see that it is possible to express G a in terms of only one 
function M 3 % k (2.29). 

The /3-functional part of 1 -gluon contribution Gp (3.11) has the same structure as LO 
kernel ( 2 . 12 ) 

-» 2 

Gp = - 2 Mf 4 (U 10 tU 02 t 3 4t + U 4 t 0 U 12 t3ot - (0 o 1 = 0 o 4)) 
no ?T 0 
-» 2 

+ 2 ^ 2 A4 (UiotU 0 2t 3 4t + U 2 t 0 U 10 t34t - (0 -) 1 = 0 O 2)) 

no ^ 20 " 

2 

— n- 9 —» 2^24 (Ul 0 tU 02 t34t + U 3 0 tU 0 4 t 12 t - (0 O 4 = 0 O 2)) 

2 ? 20^40 

2 

— 9 - 9 —* 2 ^13 (U4toUi2t30t + U 2 toU 34 t 10 t - (0 o 1 = 0 o 3)) + (1 o 3, 2 o 4). 

2 ? 10 ^30 

(3.22) 


Here AT 3 is defined in (2.33). The 1-gluon term without beta function reads 


G = Gi + G 0 . 


(3.23) 


In G one can pick the terms independent of r*o and integrate them out if they are convergent. 
We call these terms Gq. In fact the choice of Gq is not unique. We have 


Go=^(U4t 1 U2t3-U4 t3 U2t 1 ) 


x In 

x In 

x 


r i 4 


ri 0 2 r 40 2 


+ 


r-23 


r 13' 


n >4 


r20 2 ?Ao 2 


rio 2 r30 2 


no 2 Go 2 


no) 
n 2 s 
no 2 
n 4 2 
n2 2 


In 

In 


no 2 no 2 


r 20 

n 2 2 

Go 2 

?q 4 2 

F 14 2 


+ 


H3 


no 2 ^3o 2 


+ 


T 2 4 


r34 


n2 


no 2 no 2 


nio 2 no 2 


r 2 o 2 no 2 


+ In 


?’ 20 ' 


In 


no 


r w 2 no 2 


f24 J \'f24 
|(1h3,2o 4). 


+ In 


no 

?43 2 


In 


nso 

n 3 2 


(3.24) 
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It has zero dipole limit. 


G = 


r 12 ' 


no 2 no 2 


In 


no 

ri2 S 


In 


— 2 ^ (2A r c U 2 t34tl — U 0 f 1 U 2 t34to — U 2t0 U4t 10 t3) 
?12 / l 


+ (U 2 tiot34tO “ U 2 t 0 U 4 t3Uoti — (0 —)• 1))} 


+ 


ri4 


no 2 no 2 


In 


no; 

nT 


In 


no 

n4 2 


in - ( 2 -^cU 2 t 34 ti — U 0 t 1 U 2 t 34 to — U 4 4 0 U 2 43 otl) 


+ (U 2 t3oti4to — U 4 toU 2 t 3 Uot 1 — (On 1))} 


1 

+ 2 


r 13 


n 0 2 r 30 2 


In 


no 

n3 s 


In 


no 

n 3 2 


+ 


n4 


no 2 no 2 


In 


no 

n 4 2 


In 


no 

nr 


x {(U 4 toU 2 t 1 + 1144411210)11013 — U 2 t 04 t40t 3 — U 2 t04t 30 t4 -(On 3)} 


+ {U 2 toU 4 t4U 0 t 3 — U 2 t o U0t4U 3 4t + U 2 t4gt 34 to — U Q f 32 ^ 0444 } 


1 


f23 


n 2 


2 f 20 2 V no 


In 


no~ 

n 3 5 


In 


no 

n3 s 


+ 


1 


n4 


2 n 0 2 V no 


n2 

no 5 


In 


no 

n4 2 


In 


no 

n4 2 


x {U 2 t3U 4 40U O 44 — U 24 OU 044 U 344 + U 2 4iot34tO — ^ 44440243 } 


+ {U 4 toU 2 44U 0 4 3 - A^ c U 2 toU 4 44 0 43 + U 2 t 34 ti - U 2 4o4t30tl} 

i ( n 3 ^ _ nT^\ , (, ( no^\ i f n^ 

x 2r*3o 2 Vno 2 no 2 / n \n2 2 / n Vn2 2 / 2 r 40 2 Vno 2 

x {U 4 t 0 U 2 44U 0 4 3 — Af c U 0 44U 2 t344o + U 2 t 3 4ti — U 2 4o4t30tl} 



In 


no 2 

n2 2 


+ {U 2 toU 4 44U 0 4 3 - A^U 4 40U4 2 4304 + U 2 t 3 4ti - U 2 4o4tlOt 3 } 

l ( n4 2 n 3 2 \ /n 0 2 \ /n 0 2 \ l_ /n 2 2 rW 2 / . /n 0 2 \ /no 2 

X 2r*3 0 2 \n 0 2 no 2 / n Vn 4 2 / n Vn 4 2 / 2r 20 2 Vno 2 n 0 2 / n Vn 4 2 / n \fi 4 2 

x {U 2 toU 4 44U 0 4 3 - Af c U O 44U02t344 + U 2 t 3 4ti — U 2 4o4tiot 3 } + (1 n 3,2 n 4). (3.25) 

All the integrals with the functions G s , G a , Gp and G are convergent. It is clear from the 
explicit expressions for Gp and G. For G s and G a one can see it recalling that has 
unintegrable singularity at ro = n and M k has unintegrable singularity at n = n = fy 
In all expressions in this section these singularities cancel. 

3.2 Double dipole 

The symmetric contribution reads 

G s = G s 4 + n/Gg + G s2 , (3.26) 


12 - 



G S 1 — ({U 0 t 1 2 t 5 U 4 t 35 to - U 0 t 5 U 2 ti 5 t 34 to — (5 —>■ 0)} + (5 0)) (1/14 — L 13 + L 2 3 - L 2 a) 

+(1 o 3,2 o 4), (3.27) 

U 2 tiU 0t5 , 


r = 1 ( f U/tt3 

q 2 ^ N r 


( U 0tl2t5 + U 0 t 52 tl - 


N r . 


') _ U 0 ti 2 t54t3 - ( 5 0)} + (5a 0)) 


x (-^14 — -^13 + -^23 — -^ 24 ) + (1^3,2f>4), 


(3.28) 


G S 2 —^(U 0 t 54 t 3 U 2 t 05 t 1 — U 0 t 5 U 2 t 15 t 34 to + (5 <-)• 0))(M 4 + M 3 — M, — M 2 + (5-H-0)) 


+ 2^4toUot35tl2t5 + U4t 0 U 0 t 52 t 15 t3 - U 0 t 3 U 2 t 15 t 0 4t5 — U 0 t 3 U 2 t 5 4to5tl + (5 O 0)) 


x (M 4 13 + M 3 14 - Mf - M 3 24 + (5«0)) + (1h3,2« 4). 


(3.29) 


Here L, L q , and M qk are introduced in (2.22), (2.23), and (2.29). The antisymmetric 
contribution reads 


G a — 2 (Uot54tgU 2 t05tl — U 0 t 52 t 1 U 4 t05t3 + U 5 toUoti 2 t5 4 t3 — U Q f 5 U 2 t 15 t3 4 to — (5 <->■ 0 )) 
x (Mi 1 - M 3 11 + M 3 12 - M 4 12 + M 3 21 - M 4 21 - M 3 22 + M| 2 ) 

+ 2 (G5t3U 0 ti2t04t5 + U 5 t3U 0 t 5 4t 02 t 1 - U 4 t 5 U 0 ti2t05t3 — U 4 t 5 U 0 t 35 t 02 ti — (5 -H- 0)) 


x (Mf + Ml 4 - Ml 3 - Ml 4 - Ml 1 + M 4 32 - Ml 1 + M 3 42 ) + (1 h3,2h 4). 
The /3-functional contribution has the form 

G/J = 


fl 3 2 M tt _ 


mP _lii 

2^„„2 iW 23 


r Mu + J^^mP 


J10 r 30 * r 20 *r 30 

x (U 2 ti 4 t 3 + U 2t34tl — U 2 t 10 t34t0 — U 2 t 0 4t30ti), 


o 9^14 T 9 —# o lv± 2A 

r io z Uo z r 2 o z r 4 o z 


(3.30) 


(3.31) 


where M ' 3 is introduced in (2.33). The remaining contribution reads 

G = Gi + G 0 . 


(3.32) 


Gq—^( 2U 2 t 1 U 4 t 3 — A^ c U 2 t 14 t 3 — 3V c U 2 t 34 t 1 ) 


+ 


r 14 


2fi 3 2 _ny 

roi 2 ro3 2 ?% 2 r 0 2 2 ’ H)i 2 r 0 4 2 


xlnlTibnfhl) 

V r 03 2 

2 

r 13 


+ 


+ 


r or 

zf 2 
712 


+ 


H2 


+ 


r 14 


2r 24 2 

H)2 2 ?h4 2 roi 2 H)2 2 ' H)i 2 HJ4 2 




47^14 2 


7-24 


+ 


7*34 


roi 2 T 02 2 r 0 i 2 r 0 3 2 ?% 2 f 04 2 r 0 2 2 r 04 2 r 0 3 2 ?i)4 2 


7’02 
2 


nn 


ln(^ in^ 14 


7-04 

2 


ri 3 2 _H 4 2 _ t ^ 2 r 24 2 \ / n ^ 2 \ (ry^ 1 

'doi 2 To 3 2 r 0 i 2 r 0 4 2 r 0 2 2 ro3 2 ro2 2 r 0 4 2 J V^oi 2 / U \n) 2 2 


T04 

+ (1«3,2^ 4). 

(3.33) 


Gi —-(U 2 t 0 U 0 t 14 t 3 + U 2t0 U 0 t3 4tl - U 0 t 1 U 2 t 04 t3 — U 0tl U 2 t 34 t 0 ) 


In | In 


ri3 


r 12 


7-23 


Hjl 2 nj3 2 r 0 i 2 r 0 2 2 74j2 2 'do.3 2 


+ 


7-23 


T’12 


7-13 


7-02 ^ 
2 


T’03" 


ro2 2 ro3 2 r 0 i 2 ro2 2 7% 2 r 03 2 


13 - 




1 ( ff 13 2 \ 1 / ^13 2 \ , ( H4 2 \ , ( n 4 2 \ ( r\2 ri4 2 ^24 2 

11 \nji 2 J \n)3 2 J \A)i 2 J n V^04 2 ) \roi 2 r 02 2 r 0 i 2 r 04 2 r 0 2 2 r 04 2 

f n 2 2 ru 2 r 24 2 V /r 2 4 2 \, /V 2 4 2 Y 
V ro,Vo 2 2 r 01 2 r 0 4 2 + r 0 2 2 r 04 2 J “ l r 02 2 J * Uk 2 / 


+^(2U 2 tiU 4 t 3 — -/V c U 0 ti2to 4 t3 ~ -^cU 0 t 34 to2ti) In (^ 4 ] l n (- 4 ^ 
2 V 7 oi J \ r 02 

x (44 - - ^44] + (1 o 3,2 ^ 4). 

V r oi r 04 2 r 0 i 2 r 03 - r 0 2 -r 03 2 r 0 2 2 r 0 4 2 / 


(3.34) 


As for the quadrupole, it is straightforward to check that none of the functions G s , G a , 
G/ 3 , G has unintegrable singularities. 

4 Quasi-conformal evolution equation for composite operators 

To construct composite conformal operators we use the prescription [9] (see also [21]) 


O con f = Q + 


2 dr] 


where a is an arbitrary constant. The conformal dipole reads [9] 


= U 2tl + PP /In (^4] (U 2t0 U 0tl - iV c U 2tl ). (4.2) 

12 4vr Z J r w z r 20 ~ \rio 2 r 20 2 y 


The evolution equation for this operator [9] is quasi-conformal 


¥ = Y S if °S^ (* + !<) (u ^ u »" - N ^ )am ‘ 

a 2 r 

+ J 4?’orf?’5{((U 0 f5 2 t 0 5tl - UotlU2t5U5to ~ (0 5)) + (0 o 5)) 

+Lf 2 (Uot5U2toU 5 ti -(Oo 5)) 

-2nfL\ 2 {tr{t a Uit b U\)tr(t a U b t b {Ul - U\)) + (5 o 0))}, (4.3) 


where Aff 2 is defined in (2.33); Lfj = L c (fi,fj ) and Lfj = L c (fi , fj ) were introduced 
this form in [11] 


T C =T fh 2 2 , f r 02 2 ri5 2 \ n 2 2 / r 0 i 2 r 25 2 \ 

12 12 4r 0 i 2 r 0 5 2 r 25 2 n Vro5 2 ?q 2 2 y 4r 02 2 r 0 5 2 ri5 2 n V r 05 2 ri 2 2 y ’ 


4? 2,pr 2 
A)2 ? 15 


fC _ f _ ' 12 . / '02 '15 \_ ' 12 . / '01 '25 

12 12 4r 0 i 2 r 0 5 2 r 25 2 Vro5 2 n2 2 y 4r 02 2 r 05 2 ri5 2 n Uos 2 ^ 2 


^ 2^ 2 
All r 25 
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For the conformal quadrupole operator using (2.12) we have 


U 12?34t - U 12t34t + Jdfo 


/ ^14 2 , f aru 2 

X l no 2 r 40 2 n V^io 2 ^40 2 
ri2 2 , / ari2 2 \ 

no 2 r 2 o 2 n vio 2 ^ 20 2 / 
r 2 4 2 1 / Q^24 2 ^ 

2? ? 20 2 ^40 2 11 V^0 2 r 40 2 y 
n 3 2 , / ari3 2 ^ 

2rio 2 r3o 2 n \rio 2 r 3 o 2 y 
+ (1 ** 3,2 o4)}. 


^ (Ui 0 tU 02 t34t + U 4 f 0 U 12 t30t — (0^1)) 
(Ui 0 tU 0 2t34t + U 2t0 U 10 t34t - (0 ->• 1)) 

) (UiotU 02 t 34 t + U 30 tU 04 t 12 t — (0 —» 4)) 

) (U 4 t 0 U 12 t3ot + U 2 t 0 U 34 t 10 t - (0^1)) 


(4.6) 


The conformal double dipole operator reads 


(Ui 2 tU 34t ) co,i/ —U 12 tU 34 t + ^2 J ^o(U 2 t 14 t 3 + U 2 t 34 t! - U 2 t 10 t 34 t 0 - U 2 t 04 t 30 ti) 


r i3 


no 2 ^30 2 


In 


ru 


no 2 ?To 2 


In 


aru 

no 2 r 3 o 2 

an 4 2 

no 2 ^40 2 


V23 


r20 2 r 30 2 


In 


+ 


r 2 4' 


r20 2 r4o 2 


In 


ar 23 

r 20 2 r30 2 

ar 24 2 

r20 2 ^40 2 


+U 4t3 (U“f—U 2tl ) + U 2tl (U^-U 4t3 ). 


conf 


(4.7) 


The evolution equations for the conformal quadrupole and double dipole operators in the 
conformal basis have the general form 


fConf ,■ ( ^ 2 

= 4^ W^( U i° tU ° 2t34t + U4t ° Ul2t30t " 1 )) CO ' l/ 

+ -.^o- 2 (Ul 0 tU 02 t 34 t + U 2 toUiot34t — (0 -> l)) con f 
HO ^20 

— 2 - 2 (Ul0tU 02 t34t + U 30 tU 04 t 12 t - (0 —» A)) con - f 

2 r 2 o r 40 

— 0 -^ 2 - 2 ( U 4toUi 2 t 30 t +U 2 t 0 U 34 t 10 t — (0 i)) co ' l 7 + (1 o 3, 2 o 4) | 

zrio^r3o^ J 

+ §i Jdr 0 dr 5 (G c r f +G c a onf )+^ f dr 0 (G„ + G co ^), (4.8) 
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d(Ul2fU34t) COn/ f 

\ 8tt- 


dr, 


dr 0 


+ 


r 13 


no 2 ^3o 2 


4 ^ 12 " 
rio 2 r 20 2 
^ 2 

^ 23 

r 2 0 2 ^30 2 


(U 4 t 3 U 2 toU otl - iY c U 4 t 3 U 2tl ) con/ 


ru 


no 2 ^4o 2 


+ 


r 2 4' 


^20 2 ^40 2 


x (U 2 tl 4 t3 + U 2t34 tl - U 2 tl0t34t0 ~ U 2 f 04 t 30 t 1 ) con d 

+ U 34t (K nlo <S> u“f) + (193,2^4)} 

df 0 (Gp + G con f). (4.9) 


+ 


87T 4 . 


df,df S (Gr , +ar , )+^- 3 


As in the previous section, the individual NLO evolution of the dipoles here is taken out 
of the functions G 

t con f _ /• =} 2 

? 12 


( Knlo < 8 > U' 


conf v _ ^^12t 

' I2 f 7 dr, 
Therefore one can rewrite (4.9) 


a 


- / df o- 


2?r 2 J f io 2 r 20 2 


(U 2 toUoti-Af c U 2ti r n/ . (4.10) 


d(Ui2fU 3 4t) con/ 

dr. 


= U 34t - 




12 t . a s 

+ W] dr ° 


<977 


4fi2 2 

no 2 ^2o 2 


x{(U 4 t 3 U 2t0 U 0tl - fV c U 4t3 U 2tl ) co ^ - U 4t3 (U 2 t 0 u 0tl - N c U 2n ) conf } 


+ 


r 13 


r-23 


n.4 


no 2 Pio 2 


r 20 2 F 3 o 2 


ri 0 2 r 40 2 


+ 


r 24 


r 20 2 r40 2 


x (U 2 ti 4 t3 + U 2 t 34 ti — U 2 t 10 t 34 t 0 — U 2 t 04 t 3 oti) COf, ‘ 7 + (1«3,2«4)} 

+ Jdndn (OT'+G7"')+^Jdr-„ (Gg + G""'). (4.11) 

Plainly, Gg and G^ are the same as in (3.11) and (3.12). The other functions G con f will 
be given below. 

To obtain these functions one has to calculate the evolution equations for conformal 
operators (4.6, 4.7) using (2.11-2.15) and express the results in terms of conformal operators 
via (4.1). Technically, it means that one has to add to the kernels of the evolution equations 
from the previous section the corrections in the form of double integrals w.r.t. r*o and rs [9]. 
To get the conformally invariant results one has to symmetrize these corrections according 
to (3.2). Then, the terms with color operators independent of ro (or rs) can be integrated 
w.r.t. r*o (or rg) via the integrals from appendix A of [22]. Finally, the terms with color 
operators independent of both ro and r^ can be integrated with respect to both ro and rg- 
In addition to the integrals from appendix A of [22], one needs the following integral 


df 0 


f r m f 0 2 r 0 ir 0 3 


\ ivt 2 2 

V' 01'02 


^2-2 
'01'03 


Mil) 


vr 3 
= - hr 
3 


13 


12 


(4.12) 


4.1 Quadrupole 

For the symmetric contribution G™"' 7 we have 


GConf = G conf + n/ Q ? + 


(4.13) 
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where G q did not change. It is defined in (3.15) . 


QConf _ ({U 0 t 34 t 15 t 02 t 5 — Ustol^tsU^^! — (5 0)} + (5f> 0)) (lf 2 + -^32 ~ ) 

+ ({U 0 tl5t02t34t5 “ U 0 t 5 U 5 t 1 U 2 t34to — (5 — > 0)} + (5h 0)) (lf 2 + Lf - iff 
+(1o3,2o4), (4.14) 

where L c is defined in (4.4). 

G s2 / = 9 (Uot34t52t05tl — U 0 t 1 U 2 t 5 U 4 to5t3 + (5 -H- 0)) 

x (Aff 34 - Aff 43 + Aff 42 - Aff 24 + Aff 43 - Aff 34 + Aff 31 - Aff 13 ) 

+ 9 (Uot35t02t54tl _ Uot3U 2 t 5 U 4 t 15 t 0 l(5f) 0)) 
x (Aff 14 - Aff 41 + Aff 42 - Aff 24 + Aff 41 - Aff 14 + Aff 13 - Aff 31 ) 

+ 9 (Uoti5to2t34t5 — U 0 t 5 U 5 t 1 U 2 t 34 t 0 + (5 -H- 0)) (Aff 14 — Aff 41 + M f ^ - Aff 721 ) 

+ ^ (U 0 t34ti5t02t 5 - U 5 toU 2 t 5 U 0 t 34tl + (5 o 0)) (Aff 23 - A/f 32 + Aff 21 - Aff 12 ) 
+(1 3,2 **4). (4.15) 


Here A 4? jk = M c { 


Ho r'j,Vk) reads 
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=M, 
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In 


In 
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8r o i 2 fo 2 2 r*05 2 ^35 2 

zf 2zf 2 
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,2 

In 


^ 2 ^ 2zf 2 
r 05 r 13 r 3 5 


fbi 2 H)3 4 


^ 2^* 2 
TL2 7 23 


8r 0 i 2 ro 2 2 r 25 V 3 5 2 


4f 02 2 r 0 5 2 ri5 2 

In 


In 


zf 2zf 2 
7 01 ^25 

r 0 5 2 ri2 2 


—* 9 —* 9 -» 9 

7 ’orro2 7 ’23 


n 2 2 'd 2 5 2 'r;35 2 


In 


In 


fos 2 ^ 2 ^ 2 
nji 2 nj2 2 'do3 2 

rofros 2 ^ 2 ^ 2 

7 t)2 4 n5 2 A35 2 


ri2 


r 13 


8r 0 i 2 ro5 2 r35 2 


—* 9 —* 4 ~* 4 
r 01 '03 r 25 

ro2 4 fb5 2 n.3 2 Ai5 2 


+ 


4ro2 2 ro5 2 fi5 2 

zf 2 ^ 2 

7 23 r 12 

8f 0 i 2 ro2 2 r 2 5 2 r 3 5 2 


In 


In 


+_ rosW _ ln 

8f 0 1 2 ni2 2 ni5 2 2 


—* 9 —* 9 -* 4 
r 01 ^03 '25 

ro2 2 'r!j5 2 n2 2 'f'i5 2 


^ 2,-? 2 
rpi 7 25 

r 0 5 2 ri2 2 

r02 2 n2 27: 35 2 

roi 2 r 23 2 r25 2 


(4.16) 


The function Aff is conformally invariant. It does not have property (2.30). Nevertheless 
like G S 2 , G“ n/ can be rearranged to form (3.17) where instead of Aff will be Aff jfc . As 
a result G^f 4 vanishes without the shockwave. Finally, one can see that G c s on ^ can be 
formally obtained from G s via the substitution Af -4 M c , L —> fy C . 

The contribution which is the product of the antisymmetric w.r.t. 5 -H- 0 parts reads 



G conf , r-A conf 
al + Kjf a2 


+ G 


conf 

aS 


(4.17) 
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Gai^ — - (U 0 f i U 2 t 5 U 4 t 0 5t 3 + U 0 t 3 4t 52 to5tl ~ (5 O 0)) 

x (Aff 31 + Aff 13 - Aff 34 - Aff 43 - Aff 42 - Aff 24 + Aff 43 + Aff 34 - # 2134 ) 

+ 2 (Uot3U 2 t5U 4 tl 5to + U 0 t 35 t02t54tl “ (5 O 0)) 

x (Aff 13 + Aff 31 - Aff 14 - Aff 41 - Aff 42 - Aff 24 + Aff 41 + Aff 14 + f? 3241 ) 
+(1«3,2«4). (4.18) 


Here Rijki = R(fi,fj,fk,ri ) is a conformally invariant function. It reads 

o _ n2 2 , / r 02 2 ri5 2 \_ n 2 2 / rpi 2 ?^ 2 ^ 

2134 2 ?% 2 r 0 5 2 r 2 5 2 n \ro5 2 ri2 2 ) 2 r 02 2 r 05 2 n 5 2 H \' 0 ) 5 2 'r 12 2 / 

r 24 2 . / r 0 4 2 r 2 5 2 \ n3 2 ] / r 0 3 2 ri5 2 \ 

2 f 0 2 2 ro 5 2 r 45 2 n \fo5 2 r 2 4 2 ) 2 f 0 i 2 r 0 5 2 f 35 2 Vn)5 2 n3 2 / 


(4.19) 


Q C a2 f — 2 (Uot34tl5t02t5 ~ (5 O 0)) (f?231 ~ -^123 

+2Lf 3 - Aff 23 - Aff 32 + Aff 21 + Aff 12 + Aff 11 + Aff 22 - Aff 22 - Aff 33 ) 

(U 0 t 15t02t34t5 _ (5 O 0)) (i?i24 — -^142 

+2Lf 2 - Aff 22 - Aff 11 + Aff 44 + Aff 11 + Aff 14 + Aff 41 - Aff 21 - Aff 12 ) 
+(lo3,2o4). (4.20) 


Here the function L c is defined in (4.5) and Rijk = R(^i,rj,fk) is another conformally 
invariant function. It reads 


-R 123 


ffi3 2 ] f H)3 2 n 5 2 \ 

2foi 2 f 0 5 2 r 35 2 n \H)5 2 rL3 2 / 


+ 


zi 2 

7 23 

2f 0 2 2 fo5 2 r 35 2 


In 


/ rp3 2 r25 2 \ 
\r 0 5 2 r23 2 J 


(4.21) 


In fact, there is freedom in the definition of the functions Aff jfc , Rijk and Rijki since one 
can redistribute terms between them. For example, one can try to redefine Aff so that 
to make the functions R zero. 

The remaining antisymmetric contribution reads 

GZ nf =l (Uot5U 5tl U 2 t 34 to -(5o 0)) 

x (2lf 2 + 2Zf 4 - Aff 44 - Aff 22 + Aff 14 + Aff 41 + Aff 21 + Aff 12 + R 241 ) 

+ 2 ( U 5t0 U 2t5 U 0t34tl - (5 o 0)) 

x (2L 21 + 2Lg - Aff 33 - Aff 11 + Aff 32 + Aff 23 + Aff 21 + Aff 12 + R 132 ) 
+(lo3,2o4). (4.22) 


The integrated w.r.t. F 5 part of the kernel has the form 

QConf = G conf + QConf (4.23) 


18 - 



Here 


G conf 

1 


(Uot3U 2tl U 4 t 0 — U 0 t 12 t04t3 


U 2t3 4tl(^ c 2 -1)) 


X 

+ 

+ 


r 14 


roi 2 ro 4 2 


r-23 


r 0 3 2 r 0 2 2 


r 34 


r 0 4 2 ro3 2 


In 2 


In 2 


In 


2 ( r 0 2 2 n 4 

\ro4 2 ri2 

rpi 2 r 2 3 2 \ 

r 03 2 ri2 2 y 

H)l 2 r34 2 \ 
Hj4 2 n3 2 / 



2 y H)3 2 ri4 
\ro4 2 ri3 
ro3 2 f 2 4 2 \ 

n)4 2 r 23 2 J . 

ro2 2 r 3 4 2 \ 

rb3 2 r 2 4 2 J 



_ n 3 2 2 / \ 

?% 2 H)3 2 n VH)3 2 ri2 2 y 
^24 2 , 2 H)l 2 ^24 2 \ 

tq 4 2 Hj 2 2 n \ro4 2 ri2 2 / 


+(UotiU 2 t 3 U4t 0 — U 0 t 14 t 02 t 3 — U 2 t 34 tl (iV c 2 — 1)) 

2 


1 

X 4 


^34' 


H)3 2 H)4 2 


In 2 


i? 2zf 2 
74)2 ^34 

m 2 ^ 2 


-In 2 


^ 2zf 2 

Hn ^34 
nj4 2 n 3 2 


?’24' 


ro2 2 ro4 2 


In 2 


zf 2zf 2 

74)3 7~24 
ro4 2 r 23 2 


+ 


+ 


r 12 


H)2 2 roi 2 

ri4 2 

r 04 2 roi 2 


In 2 


In 2 


zf 2zf 2 

74)1 ^23 
n)3 2 n 2 2 

zf 2zf 2 
7 03 ^ 14 

nj4 2 n 3 2 


-In 2 


+ In 2 


zf 2 zf 2 

74)1 7 24 

f 04 2 ri 2 2 

2^ 2 
? 01 ^24 

ro2 2 r 14 2 


r 13 


H)3 2 Hji 2 


In 2 


^ 2 4? 2 

74)1 7 23 

r 0 2 2 ri3 2 


+ (1 3,2 o 4). 


(4.24) 


G 2 ^ — ~^ _ (U2toU 0 t3 4 ti — ^cU 2 t3 4 ti) 


7*13 


roiVos 2 


In 2 


zf 2zf 2 

74)2 7-13 


na 2 ^ 2 


r-23 


7 ? 02 2 n)3 2 


In 2 


7"oi 2 r23 2 

?4)3 2 ri2 2 


T’12 


7% 2 H)2 2 


iV, 


+ “^(Uot3U 2 t04tl - iV c U 2t34tl ) 


7~24 


In 2 

2 


^ 2,-? 2 
74)1 7 23 


Hj3 2 ri2 2 


HJ2 2 H)4 2 


In 2 


zf 2zf 2 

rp3 7-24 

r 0 4 2 r 23 2 


7-23 


r 0 2 2 f03 2 


In 2 


zf 2zf 2 

7 02 7-34 
7 ? 04 2 7 ; 23 2 


7’34 




In 2 


zf 2 zf 2 

74)2 7 34 
7 ? 04 2 7 ; 23 2 


+ In 2 


+ In 2 


zf 2zf 2 
7 01 7-23 

7^)2 2 74 3 2 


2^ 2 
7 02 7-34 

?4)3 2 r 24 2 


+ (1 3,2 4). (4.25) 


It was straightforwardly checked that all the integrals of G c s on ^, Ga° n f and G con f do not 
have unintegrable singularities. 


4.2 Double dipole 

For symmetric contribution Gf"^ we have 

Gf n/ = G c s ° nf + n f G q + G“ n/ , (4.26) 

where G q (3.15) did not change. 

Qcon/ _ ({U 0 ti 2 t 5 U 4 t 35 t 0 _ U 0 t 5 U 2 t 15 t 34 t 0 — (5 —■> 0)} + (5f> 0)) (if 4 — if s + L 23 — if) 
+(1 <-)• 3, 2 -H- 4), (4.27) 

G c s ° nf = (M 4 C13 - M 4 C31 + Mf 14 - Mf 41 - M 4 C23 + Mf 32 - Mf 24 + Mf 42 ) 

X 2 (^4toUot35tl2t5 + U 4 toU 0 t52ti5t3 ~ U 0 t 3 U 2 t 15 t 0 4t5 — + (5f> 0)) 

+(M 4 C12 - Mf 21 + Mf 21 - Mf 12 - Mf 34 + Mf 43 - Mf 43 + Mf 34 ) 

x 2 (^0754+3^2105+1 _ U 0 t 5 U 2 t 15 t 34 t 0 + (5 0)) + (1 3,2 n 4). (4.28) 
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The antisymmetric contribution reads 


Ga — 2(^5t3Uotl2t04t5 + U 5 t 3 U 0 t 5 4to2t 1 “ U 4 t5U 0 t 12 to5t3 — U 4 t 5 U 0 t3 5 f 02 ti -(5f> 0)) 

X (Mf + M 3 24 - M | 3 - M 3 14 - Mf + Mf - M 3 41 + M 3 42 - i?341 + i?342) 

+ 2 ^° t54t3 ^ 2t05tl ~~ ^0t52tlU 4 t05t3 + U 5 t 0 U 0 t 12 t 54 t 3 — U 0 t 5 U 2 t 1 5t 34 to — (5 o 0)) 

x (Mf 11 - M^ 11 + A/f 12 - Mf 12 + M 3 C21 - Mf 21 - M 3 C22 + Mf 22 ) + (1 o 3,2 o 4). 

(4.29) 


The contribution with 1 gluon crossing the shockwave has the form 


G conf — ^(U 2 t 0 U 0 ti 4 t 3 + U 2 t 0 U 0 t 34 t 1 U 0tl U 2 t 04 t 3 U 0tl U 2 t 34 t 0 ) 


r 12 


rm 2 rQ 2 2 


In 2 


z* 2zz 2 

74)2 ^14 


r 0 4 2 ri2 2 


-In 2 


zf 2zf 2 

7 02 ^ 13 


'd)3 2 ri2 2 


+ 


ri3 


r23 


roi 2 ro3 2 


ro2 2 r 0 3 2 


In 2 


roi 2 7"23 2 

ro2 2 7h3 2 


-In 2 


r 14 " 

rolVo4 2 


i? 2^ 2 

74)1 7-23 

ny.i 2 r\ 2 2 


+ 


r24 


ro2 2 ro4 2 


+ In 2 
In 2 


zf 2z* 2 

74 ) l t - 24 
r 0 4 2 ri2 2 
rpi 2 r 2 2 
F 02 M 4 2 


+ 4 (N C XJ 0+12+04+3 + -^cU 0 t 34 t02+l “ 211213X1341) 


+’24" 


ro2 2 r 0 4 2 


In 2 


+ 


ri3 


r 0 i 2 f 0 3 2 


In 2 


zf 2zf 2 

74)2 r\3 

r 0 3 2 ri2 2 


ru 


fo i 2 r 04 2 


In 2 


zf 2zt 2 

r 02 ri 4 

r 0 4 2 ri2 2 


+ In 2 


zf 2 rf 2 
74)2 +34 

nj3 2 '+ 7 24 2 

2,^* 2 
74)1 7 34 

r 0 3 2 744 2 


+ (1 3,2 ^ 4). 


(4.30) 


As for the quadrupole, it was straightforwardly checked that all the integrals of G“ n , 
Ga , and G con ^ do not have unintegrable singularities. 


5 Checks 

There are two checks which can be done for the results of this paper. The evolution 
equations for the quadrupole and double dipole operators can be obtained from the NLO 
JIMWLK hamiltonian [16] and the general evolution equations from [17]. 

In this paper the following two checks were done. First, quadrupole kernels (3.11) and 

(4.8) go into BK ones (2.34) and (4.3) in the dipole limits 1 —)• 2, 2 —> 3, 3 —)• 4, and 4 —> 1. 
Double dipole kernels (3.12) and (4.9) also have the correct dipole limits 1 —>• 2 and 3—^4. 
In these limits they also go into the BK ones (2.34) and (4.3) times N c . This statement 
can be checked straightforwardly going to the dipole limits in explicit expressions (3.11), 

(4.8) , (3.12), and (4.9). Our kernels match the Balitsky-Fadin-Kuraev-Lipatov NLO kernel 
[23] in these limits. 

The second check is that in SU(3 ) our kernels respect the identity 

-B 123 = U 12 tU 34 t — XJ 12 t 34 t, (5-1) 

where B 123 is the 3-quark Wilson loop (baryon) operator defined as 

£123 = Ci • U 2 • U 3 = e^'b'sijhU^U^U^,. (5.2) 
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The evolution equations for the l.h.s. of (5.1) in the standard and quasi-conformal forms 
are given in [11], They read 


dB 


123 




ajj^_ 

87 r 2 


df 0 

3a. 


(- 6100^320 + -B 200 -B 310 — B 300 B 210 — 61?123) 


r 2 

7 12 _ on g 

47r 


2-2 
' 01'02 


InlMA 

111 1 ^2 I \ ^2 -2 -2 -2 


a. s , for, , 74o f 1 


in ' 2 ? in 


7 r r, 


21 


21 


02 


13 


r 2 
7 12 


lnl^f 




' 32 


1 2 ryf 2 (f>*2 />T 2 1 

L'l0'30 ' 30' 20 J 


(^ 100-^320 — ^200-8310) 


12 


95i23 - - [2 (-B 100 -B 320 + -B 200 -B 130 ) - -B 300 -B 120 ] ) +(1h3) + (2h3) 


7 10'20 v 


a 2 s nf 

167T 4 


df 0 dr 5 


\{UxU^U b + t/st/o^i) -U 2 -U 3 - ^B 123 tr(UjU b ) 


+ (U\U^U 2 ) ■ U;i ■ U 3 + -B 123 ~ y(i?013^002 + Bqq\Bq 2 3 - £> 012 - 6003 ) 

O 4 

+ (1^2)) +(0o 5)}L? 2 + (1 fi 3) + (2 ft 3)j 

O 

{T 12 (U 0 UJU 2 ) ■ (fWt/s) ■ U 3 

(u 0 uju 2 ) ■ (ihuju^) ■U 3 + tr ( U 0 U 5 f ) ( lhuju 2 ) • U 3 • U 5 


- W 


+T 12 


3 1 

“ ^[-^155^235 + 62556135 ~ ^355^125] + 2^ 123 


+ (M 13 - M 12 - M 23 + M 2 13 ) [(tio^t/a) • (Wo f £/i) • 65 + (T/it/o+C/s) ■ ( 6365 * 60 ) ■ 65 


+ (all 5 permutations 1 « 2 « 3)} + (0f>5) 


(5.3) 


where L, L q , M? k and M*j are introduced in (2.22), (2.23), (2.29), and (2.32).(5.3) and 
(5.4). 


dB 123 f _ a * 


(*?) 


di] 


8 tt 2 


df 0 


((B±ooB 3 20 + B 200 B 3 10 — 63006210 ) - 6 B l 23 r nf 
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r 0 2 i?o 2 4lr 


In | M W i_L | _ 
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02 

T32 2 ] 2 f ^Mp 2 
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+ (1 3) + (2 ■<-> 3) 


7 h 2 2 ] 2 / ri2 2 f 30 2 

'fin 2 '7112 2 n Vri3 2 r2o 2 


+ (all 5 permutations 1 44 2 !->• 3)^ 


a 2 rif 
167T 4 . 


df 0 df 5 


\(U 3 UjU b + 6561 * 6 ^ • 6 2 • 6 3 - l -B l 23 tr{UjU 3 ) 


+( 6160 * 62 ) ■ U 3 ■ 65 + -6123 - 7(60136002 + 60016023 - 60126003 ) 

o 4 

+ (1 H 2)^ + (0 ft 5)| + (1 <4 3) + (2 3) 
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f dr 0 dr 5 ({Zf 2 (u 0 U 5 'U 2 ) ■ (EWcM) • U 3 

+Lf 2 [( 'u q uJu 2 ) • (EfiEEo^) • £4 + tr ( U 0 U 5 f ) (EWo^) • ££3 • U 5 
3 1 

— ^[-Bi 55-££235 + -B255-Bl35 — -B355-£?125] + 2^123 

+m£ [(f/of/s^s) • (WtM) • ££5 + (eEiEEcM) • (e^EEo) • ££5 

+ (all 5 permutations 1 o2<r-3)| + (0<r- 5)^ , 
where L c is defined in (4.4), L c — in (4.5), and 


M ' 2 16r 0 2 2 ro5 2 ri5 2 ^ 

£* 23 2 ] 

16r 02 2 ro5 2 r35 2 
H3 2 , 

16r 0 3 2 £o5 2 £i5 2 


—* 2 ~* 2 ~* 4- \ —* 2 

r'01 'f'Q2 r 35 \ _ ri2 ^ 

ro3 4 £i5 2 £25 2 / 16r 0 iVo5 2 r 2 5 2 


—* 4- —* 4 —» 4 -* 2 
n)3 H)5 H2 ?" 25 

£*01 2 £02 6 £l5 2 £*35 4 
/ —* 2 ~* 2 -» 4 

/ r 0 2 r 03 ri 5 4 


£bl 4 £b3 2 £25 6 £35 2 \ £*23 2 ln / £t)2 2 £b3 2 ' ^ 

do 2 2 r 0 5 4 r 1 5 4 r 2 3 4 / 16r 03 2 ro5 2 r 25 2 VEoi 4 d 2 5 2 r 35 2 

n )2 ri 5 z r 35 z \ ^13 ln / ^02 n 5 z r 35 z 

doi 2 r 03 2 r 25 4 / 16 r 0 i 2 r 0 5 2 r 35 2 \doi 2 do 3 2 r 25 4 


doi 2 r 03 2 r 25 4 


do 3 2 ?4 2 2 / £bi 2 fb 3 2 d25 4 \ r 23 2 ?4 2 2 / £o 2 2 di 2 2 d 35 2 

8 r 0 i 2 r 02 2 ro 5 2 r 3 5 2 H vdo 2 2 do 5 2 ri 2 2 r 35 2 / 8 ?% 2 r 02 2 r 25 2 r 3 5 2 H Vr 0 i 2 r 23 2 r 25 2 

Ei 5 2 r 23 2 / 7 % 2 r 05 2 r 23 2 r 25 2 \ 

8r 0 i 2 ro5 2 r 25 2 r 3 5 2 H \ r 02 4 di 5 2 r 3 5 2 ) ' 


In order to check identity (5.1) one has to rewrite the evolution of its l.h.s. in the same 
operator basis as the r.h.s. To this end one can use the SU(S) identities 

0 =[££ 0 • Efi • U 2 tr (££ 0 ^ 5 ) tr (££ 5 t ££ 3 ) 

-tr (E/ 5 ££o t ) (U\U^U 3 + U 3 UjU x ) • ££ 0 • ££ 2 + (e/qEVeEi) • (iWo^s) ■ E £ 2 

+ (e/iE/s^o) • (E£ 5 E£o t E£3) • E £ 2 + (1 ^ 2)] - (5 ** 0), 

0 = 2 tr (E/sEV) (u 2 U^U 3 + E£ 3 E£ 5 t E£ 2 ) • U 0 ■ Ui 

+ (E£ 0 E£ 5 t E£i + E£ 1 E£ 5 t E£o) • (E£ 2 E£o t E £3 + U 3 U^U 2 ) ■ U 5 

+ (E£oE£5 t E£ 2 - EE 2 EE5 t EEo) • (E£ 3 E£o t E£i - EEiEEcM) ■ U 5 

+ (u 0 U^U 3 - E£ 3 E£ 5 t E£ 0 ) • (E£ 2 E£ 0 t E£i - E£iE£o t E£ 2 ) • E £ 5 - (5 o 0 ) 

in the antisymmetric color structures and then 

Ui • Uj • U k = ( UiU1 ) • (££,■££/) • {U k uj ) = (t//t/*) • (£// Uj) • (£//E4) (5.6) 

with l = 2 to express t / 2 in terms of I/j. After that one can expand the products of 
Levi-Civita symbols as 

I 8 ? 8{ 6 f I 


e«he‘ ,fh ' = Si S 1 ! Si 


SI si 


and see that (5.1) is satisfied. 
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6 Discussion and conclusion 


This paper presents the evolution equations for the double dipole and quadrupole operators 
in the standard (3.11), (3.12) and quasi-conformal forms (4.8), (4.9). They have correct 
dipole limits and in SU(3) obey group identity (5.1) with the corresponding evolution 
equations for the 3QWL operator obtained in [11]. This fact ensures the correctness of all 
the 3 results. To construct the composite operators, prescription (4.1) was used. It was 
proposed in [9] for the dipole and proved in [21]. Here it gave the quasi-conformal kernels for 
both double dipole and quadrupole operators, thus being checked by the specific calculation 
of the evolution of the 4-point operators. 

Unlike the dipole and the 3QWL operators, the evolution of the quadrupole and the 
double dipole ones generates several operators in the virtual part. Indeed, the virtual gluons 
do not change the color structure of a dipole or a baryon. New color structures appear in 
the evolution of these operators only when the gluons cross the shockwave. Therefore, one 
can write the virtual part of the evolution equations for them without calculation. The 
double dipole and the quadrupole, on the contrary, mix in the virtual part with each other 
and with the double dipoles and quadrupoles with the other order of the Wilson lines. 
Therefore they had to be calculated explicitly. Using the evolution equations for Wilson 
lines from [12] in this calculation, one encounters ill-defined integrals which were treated 
here so as to obtain the known result for the dipole and the 3QWL operators. Although 
such treatment gave the equations satisfying all the checks, it is important to have the 
initial expressions with the regularization of the IR singularities and to check the results 
of this paper. Such checks may be performed starting from the evolution equations found 
in [16] and [17]. 

The equations obtained in this paper may be used to derive the NLO evolution equation 
for Weizsacker-Williams gluon distribution. This work is in progress. They can also be 
important in the analysis of higher than dipole Fock components of the virtual photon in 
the diffractive processes. 
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